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Abstract. In this article we study a piecewise linear discretization schemes 
for transfer operators (Perron-Frobenius operators) associated with interval 
maps. We show how these can be used to provide rigorous pointwise approx- 
imations for invariant densities of Markov interval maps. We also derive the 
order of convergence of the approximate invariant density to the real one in 
the L^-norm. The outcome of this paper complements rigorous results on L 1 
approximations of invariant densities |14| and recent results on the formulae 
of escape rates of open dynamical systems I16| . We implement our compu- 
tations on two examples (one rigorous and one non-rigorous) to illustrate the 
feasibility and efficiency of our schemes. 



1. Introduction 

Although this article is about the approximation of invariant densities for interval 
maps, it is intimately related to what are commonly termed open dynamical systems 
or maps with 'holes' [BJ. Open dynamical systems have become a very active area of 
research. In part, this is due to their connection to metastable dynamical systems 
[TU1 [TTj and their applications in earth and ocean sciences [5J [TB] . Corresponding 
to invariant measures for closed dynamics, in open dynamical systems, long-term 
statistics are described by a conditionally invariant measure and its related escape 
rate, measuring the mass lost from the system per unit time [BJ. 

In their recent article [16j , Keller and Liverani obtained precise escape rate for- 
mulae for Lasota-Yorke maps with holes shrinking to a single point. These formulae 
depend, pointwise, on the invariant density of the corresponding closed system. Un- 
fortunately, explicit formulae of invariant densities for Lasota-Yorke maps are, in 
general, unavailable. Thus, to complement the result of it is natural to consider 
numerical schemes which provide rigorous and computable pointwise approxima- 
tions of invariant densities. 

In the literature, rigorous approximation results are available in the L 1 -norm |14j . 
the X°°-norm [5] and in the BV-tiottii, the space of function of bounded variation, 
[l2l [7], for example. None of these methods are well-suited to our problem. For 
example L 1 approximations cannot provide the pointwise information necessary 
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for application of the formulae of [16] (see Section 7 in 1J). Convergence schemes 
in the L°°-norm and BV-norm cannot provide a computable error bound as the 
constant in their approximation error depends on the invariant density itself, which 
is a priori unknown. 

By rigorous and computable approximation we mean the following. Assume 
we are given a transformation r (typically, a formula) and an error tolerance, for 
example A := 10~ 2 . We choose a suitable discretization scheme for the transfer 
(Frobenius-Perron) operator associated to t and we are asked to determine an 
explicit level of discretization e such that the approximate invariant density /* for 
the discretize operator at level e satisfies 

nr-/;iioo<A 

Here /* is the invariant density for r. We emphasize that we assume the continuous 
density /* is unknown throughout this calculation. By computable we mean that, 
at each step, one can determine via an algorithm, within a finite num ber o f steps, 
each quantity necessary to determine^ /* and to guarantee Inequality (1.1). 



One novelty of our approach is that we will need to consider two different dis- 
cretization schemes in order to carry out this task, both based on binned discretiza- 
tion of the state space. The results in this paper enable us to compute the number 
of bins m, with e = m , and the associated approximate density, as usual denoted 
which achieves the tolerance A, uniformly. 

The first ingredient in our analysis uses a natural piecewise linear discretization 
scheme and the abstract perturbation result of |15j . The two Banach spaces involved 
in our computation are L°° , and the space of Lipschitz continuous functions on the 
unit interval. The same Banach spaces were used in [9] to provide a computer- 
assisted estimate on the rate of decay of correlations. However, the discretization 
scheme which was used in [5] is the traditional Ulam method. Ulam's method 
does not fit in our setting since it does not preserve the regularity of the space 
of Lipschitz continuous functions. The idea of our method is to first prove an 
appropriate Lasota-Yorke inequality for transfer operators associated with Markov 
interval maps, then to construct a discretized transfer operator which preserves the 
regularity of the space of Lipschitz continuous functions and which is close, in some 
suitable norm, to the original transfer operator. Although, neither the original 
transfer operator nor its discretized counterpart is a contraction in the L^-norrrPJ 
we obtain quasi-compactness of the original transfer operator, thanks to the result 
of [13] , We use the general setting of [15] which allows the study of perturbations 
of transfer operators which are not necessarily contractions on either Banach space. 
With this machinery, we can compute an explicit upper bound on the norm of the 
resolvent, bounded away from the spectrum, of the transfer operator associated 
with the map. Once an upper bound on the norm of the resolvent is computed, we 
use a second discretization scheme, whose associated finite rank operator is Markov, 

■'"Of course, the efficiency of such an algorithm is an important issue. For our purpose, we will 
be satisfied with algorithms that can be implemented in standard mathematical software on a 
personal computer. Beyond that, we do not specifically address computational efficiency in this 
article 

2 Unlike the L 1 setting where the norm of the transfer operator is automatically \\C\\i < 1, 
here we can only show that \C\ < M, where M is typically greater than or equal to 1. 
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to compute an approximate invariant density with the pre-specified error tolerance 
A. 

The reason for using two different discretizations in our method is the following: 
The first discretization has the projection property. This property is essential in 
the proofs related to the computation of an an explicit upper bound on the norm of 
the resolvent. Moreover, it produces reasonable constant^] which are needed when 
using the perturbation result of [15] . However, this natural discretization leads to 
a non-Markovian finite rank operator. The lack of the Markov property makes the 
(theoretical) rate of convergence slow. Thus, at the next stage, we use a different 
discretization, which lacks the projection propertjQ but produces a finite rank 
operator which is Markov. With this Markov scheme we will obtain a computable 
rate of convergence which is of order to -1 In to. 

In Section [2] we set up our notation and assumptions. We also recall known 
results on Markov interval maps which are needed in the sequel. In Section [3] we 
provide a Lasota-Yorke inequality on the space of Lipschitz continuous functions. 
In Section[4]we present our two discretization schemes and prove results about their 
regularity properties when acting on the space of Lipschitz continuous functions. 
In Section [5] we present the perturbation result of |15j as a sequence of steps which 
are necessary for rigorous computations and in Section [6] we apply the perturbation 
result to our problem. The main challenge of this paper lies in this section where 
we design an algorithm which enables one to rigorously compute an upper bound 
on the norm of the resolvent of the continuous transfer operator. The resolvent 
estimate can then be used to obtain a rate of convergence CrnT 1 him in Section 
[7J where C is a computable constant. In Section [8] we implement the algorithm of 
Section [6] on a specific Markov interval map. Numerical results are reported for all 
critical constants. In particular, we compute the number of bins (to = 7 x 10 6 ) 
that guarantee approximation of /* by f* n within tolerance A = 10~ 2 . Section [9] 
contains a discussion of an alternate projection based (non Markov) scheme and 
in general, on efficiency of both piecewise linear discretization schemes for uniform 
approximation. 

2. Assumptions and Notation 

2.1. Function spaces. Let (1, 03, A) denote the measure space where / := [0, 1], 03 
is Borel u-algebra and A is Lebesuge measure on /. Let L p = L' P (I) = L P (/,03,A) 
for 1 < p < 00. Let C Lip (7) denote the space of Lipschitz continuous functions on 
/. We equip C Llp (I) with the norm 

IHI=Lip(-) + H, 

where Lip(-) is the Lipschitz constant of a function and | • | is its L°° norm. The 
Arzela-Ascoli theorem implies that the unit ball of || • || is | • |-compac1n 



3 This is very important from computational point of view. In particular, smaller constants 
means that less time will be spent on the computer to produce the desired computation. 
4 Thus, we could not use the Markov discretization right from the beginning. 
^See [8| for this and other basic properties about the Banach space C Llp (7). 
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2.2. Markov maps of the unit interval. Let r : I — > / be a measurable trans- 
formation. We assume that there exists a partition V of /, V — {ci}f =0 such that 

(1) for each i = 1, . . . , q, Ti := Ti ( Ci _ 1 , Ci ) is monotone, C 1+hlp and it extends to 
a C 1 function on [c,_i,Cj]; 

(2) (t^~ )' : [r(ci_i), r(cj)] — > [ci~i,Ci] is a Lipschitz continuous function for 
each i. 

(3) there exits a number ao such that rpj < cto < 1; 

(4) for each i, j = 1, . . . , q, if r((ci_i,Cj)) n (cj-x,cj) ^ 0, then r((ci_i,Ci)) 3 
(Cj_i, Cj); 

(5) there exists a fco > 1 such that r k ° is piecewise onto. 

Let C T : L 1 — > L 1 denote the transfer operator (Perron-Frobenius) [2] associated to 

y— r L x 

2.3. Results implied by the assumptions on r. Assumptions (l)-(5) imply (see 
[lj) that C T : Cfo?{I) —> the space of locally Lipschitz functions (relative 
to the Markov partition V) and C T k : C Lip (I) -t C Lip (I) whenever k > k . Fur- 
thermore, r admits a unique absolutely continuous probability measure dv — f*d\ 
equivalent to Lebesgue, with /* G C Llp (/) and uniformly bounded away from zero 
([3] Theorem 0.1). Further properties of the system {t,u) are as follows: 

(1) The system is Bernoulli (;4s Theorem 0.3) hence t is mixing with respect 
to v. 

(2) The eigenvalue 1 is a simple eigenvalue for the operator C T acting on cither 
of the Banach spaces L 2 (I, v) or C, 'i'(J). 

(3) 1 is the only eigenvalue of unit modulus for the operator C T acting on cither 
Banach space. 

A useful regularity estimate playing a critical role in the proof of Theorem 0.1 in 
[4] is the following: 

(2.1) B := sup Lip(£?l) < oo 



3. Lasota-Yorke inequality on C Uf (I) 
Lemma 3.1. For all n > 1 

\£»\ < M, 

where M = B + 1 . 

Proof. Since each function C n l is nonnegative on / with unit L 1 — norm, it follows 



that |£"1| < 1 + Lip(£ n l). The result now follows from (2.1 1. □ 



We now fix a fc > fco such that a = Mo>q < 1. Let T := r k and C — C k be the 
transfer operator associated to T. We also set To = {ai}\ =0 to be the partition of 
/ which consists of the intervals of monotonicity of T. 

Remark 3.2. By the mixing property of (r, v) r and T have the same unique invari- 
ant density /* . While our goal is to provide a rigorous approximation of the density 
/* for t in the L°° norm, we will achieve this goal by working with T instead of r. 
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Lemma 3.3. We have the following Lasota-Yorke-type inequality for C : C Llp (I) — > 
C Llp (I) and 

Lip{Cf)<aLip{f)+B x \f\, 
where B 1 := Lip(£ k T l) (= Lip(£l)). 

Proof. Let x,y e I. For / e C Lip (J) we have 

< ' f( T r 1 x )~f(Tr 1 y) ' /prV) /(ITSO 

* Li p(/) El ^feSr^ l + l/l E ICO'OO - PTWl 

i= 1 I V i /I i=l 

- i 

^ «o E M Li p(/)l a; ~y\ + Li p( £ i)l/ Ik - wl 



<aLip(/)|.T-y|+B 1 |/||.x-y| 



□ 



Corollary 3.4. For f 6 C Llp (I) and for all n > 1, we have 

||£"/||<a"||/|| + M(l + -^)|/|. 

1 — a 

Proof. The proof follows from Lemmas |3.3| and |3.1| □ 



Remark 3.5. The inequality in Corollary |3.4[ combined with the mixing property 
of T gives the following spectral picture for £ acting on C Llp (I): 

(1) The essential spectral radius^] of C is p css < a < 1 

(2) The spectrum of C outside of the disc {p < a} consists of a simple eigenvalue 
at p — 1 corresponding to the unique invariant density for T (or r). 

This picture will play a critical role in what follows. 

3.1. Estimates on the relevant constants. The key expression from the previ- 
ous section is the inequality in Corollary |3.4| Reviewing the constants that need to 
be estimated, a and B\ are easily computed directly from the definition of t, but 
the constant M, which is defined in terms of Bq may not be so straightforward. 
One favourable situation is when a variational Lasota-Yorke inequality is available 
for r. In that case we can estimate M directly as follows. We assume C T satisfies 

(3.1) V£»f< A a n Vf + B\\f\\, n>l 

where V(-) denotes the usual notion of variation for L 1 functions and || • || the L 1 — 
norm. Then we simply observe 

< |£™i| < v{£»i) + \\c*i\\ < A a n vi + 511111 + 1 = 5 + 1 

so we can set M := B + 1. 

''See, for example Hennion 1131 for an elegant proof of this fact. 
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4. TWO PIECEWISE LINEAR DISCRETIZATION SCHEMES FOR £ 

4.1. Projection-Based Discretization. Let r\ — {bi}™^ be a partition of I into 
intervals of size at most e; i.e., mesh(?7) < e. Since uniform partitions are the 
first choice for numerical work, we will assume for the rest of this article that 
bi — bi-i = — . Everything we do could be modified for non-uniform partitions with 
only minor notational changes. Let 

V% = Xpn-iA]'^ = t / Pid^ 1 = -, \ — = m and <t>i(x) = m [ <Pid\- 

JO Oi ~ JO 

For / : I -> K define 

m 

n m f=f(b ) + J2Uf(h)-f(bi- i )} 

i=l 

m— 1 

= /(6 )(1 - 0i) + f{b m )4> m + ~ 

Thus, II TO / is a piecewise linear function with respect to rj; moreover, Tl m f(bi) = 
f{bi), i = 0, . . . , m. We shall write fl m in a more compact form 

m 
i=0 

where the basis ipi is a set of hat functions over r\: 

(4.1) ip := (1 - <f>\) ,i>m ■= 4>m and for i = 1, . . . , m - 1 , ipi := (fa - fa+x). 
Lemma 4.1. The operator H m is a projection; i.e., 

TT 2 — TT 

Proof. Observe that 

1 if i = j 



For / G C Up (I), we have 



m / rn 



n m (n m /) = E f( b i)^( b j) ^ = E /foWi = n " 

j=0 Vi=0 / j=0 



□ 



Remark 4.2. The projection property of II TO is going to be essential for results in 
Section [6] (specifically, in the proof of Lemma 6.3). 

Lemma 4.3. For f G C Llp (I) we have 

(1) Lip(Tl m f) < Lip(f); 

(2) |H m /| < |/|; 

(3) |II m / - /| < 2 £ £*p(/); 

(4) vn m / < vf. 

Proof. For / G C Lip {I). We have 

Lip(n m /) = max = max < Lip(/), 

i bi — % bi — bi-i 
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and 

m m 

|iW(*)l = |^/(6i)^(x)| < = l/l- 

2 = 1 = 

To prove the third statement, let x £ [b i _ 1 ,b i \. Then 
|n„J(x) - f(x)\ = |n m /(z) - n m f(bi) + f(h) - f(x)\ 

< \n m \\f(x) - f(bi)\ + - f(x)\ < 2\f(bi) ~ f(x)\ < 2eUp(f). 

Finally, we have 

m— 1 

VU m f=Y / \f(bi+i)-f(bi)\<Vf. 

□ 

We now define a discretized version of £ by: 

C — TT rn 

Observe that C m is a finite rank operator whose range is contained in the space of 
continuous, piecewise linear functions with respect to rj. On this space, and with 
respect to the basis {ipi}, the action of C m becomes matrix multiplication^] by an 
(m + 1) x (to + 1) matrix whose (ij) th entry is given by 

m i: j := £ipi(bj). 

By definition, for all i, the functions tpi are non-negative. Consequently, the matrix 
corresponding to C m is a non-negative matrix. However, the matrix of C m is not 
necessarily stochastic. 

4.2. Markov Discretization. Another natural choice for constructing piecewise 
linear approximations was used by Ding and Li 7 . In this case the goal is slightly 
different, namely, to preserve the Markov structure of the tran sfer operator C. 

Let ipi, i = 0, 1, 2, . . . to be the hat functions defined in (4.1 ). For / € L 1 , we set 
h ■= [h-i,h] and 

fi := to J fdx, 2 = 1,2,... to, 
the average of / over the associated partition cell. For / G L 1 we set 



m — 1 



Qmf ■= Mo + ^2 % + 2 %+1 4>i + f m lpm 
i=l 

While Q m fails to be a projection operator, it retains good stochastic propertied 
0: 

• Q m : L 1 — > L 1 is a Markov operator whose range is contained in the class 
of continuous, piecewise linear functions with respect to the bi. 

• VQ m f<Vf 

Moreover, Q m has nice regularity properties when acting of the the space C Llp (I) 

Lemma 4.4. For f e C Ltp (I), we have 
(1) \Q m f\ < 2|/|; 



7 A11 our matrices act by multiplication on the right, in keeping with the usual convention in 
the literature. 

^See [7] for proofs. 
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(2) Lip{Q m f) < lLip(f); 

(3) \Q m f-f\<^Lip(f). 

Proof. Let / E C Lip (I) and x e h, i = 2, . . . m - 1. Then, 

Qmf(x) = \ fi -^ h ^{x) + fi + * i+1 M*)\ < 21/1- 

For the special cases when x € I\ or in / m , the proof is similar. For the second 
statement we have 

Lip(Qm/) = to max - + - | = ^max|/ i+ i - £_x| 

2 z Z Z i 

= TT max l/(^+i) - /(^i-i)l) 

for some G and € Ij— i- Therefore, 

Lip(Qm/) = ^Lip(/)max|z l+ i - Zj_i | < ^Lip(/). 
z % 2 

To prove the last statement, we let x € li and observe that 
\Q m f(x) - f(x)\ = \ fi -^ fi ^{x) + fi+ 2 fi+1 M*) - f(*)\ 
= if + ^i-i{x) + f ^fUx) - 

= + + f ^fUx) (^i-i(x) + ^(*))^| 

< | - /(a?) | + | /i-i ~ /(g) | + | /i+i ~ f(x) i 

< ~[-Lip(/) + -Lip(/) + -Lip(/)] = ~Lip(/). 
2 m m m 2m 

□ 

We now define a second piecewise linear Markov discretization of C by 

]P m . Qm O 

Notice that P m is a finite-rank Markov operator whose range is contained in the 
space of continuous, piecewise linear functions with respect to rj. The matrix rep- 
resentation of P m restricted to this finite-dimensional space and with respect to the 
basis {ipi} is a (row) stochastic matrix, with entries 

Pij :=m I Ltpi > 0. 

5. Keller-Liverani's theorem 

In this section we present a version^] of the perturbation result of [15] which is 
designed for rigorous computations. All the constants involved in the statements 
below are crucial in our work. 



"in |17| another version of 1 151 was designed for rigorous computations in the framework of 
(BV, || ■ ||) and (L 1 , | ■ |), where the | ■ |-norm of both operators in the setting of |17| is smaller than 
1. 
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5.1. Notation. Let (B, || • ||) be a Banach space which is equipped with a second 
norm | • | such that 

(5.1) H<IH|. 

For any bounded linear operator P : B — > B, consider the set 

V s , r (P) = U G C : \z\ < r or dist(z, cr(P)) < 5}, 

where <r(P) is the spectrum of P with respect to (B, || ■ ||), and define 

H s , r {P) ■= su P {||(2 - P^H : z £ C \ V s , r } < oo. 

Further, we define the following operator norm: 

|||P||| = sup \Pf\. 
Il/ll<i 

5.2. Assumptions, let Pj : B — > B be two bounded linear operators, i = 1,2. 
Assume that: there is a C\, M > such that for all n E N 

(5.2) |P"|<CiM n ; 
and 3a G (0, 1), a < M, and C 2 , C 3 > such that 

(5.3) ||if/||<C 2 a"||/||+C 3 M"|/| VneN V/el,i = l,2; 
moreover, if 

(5.4) \z\ > a, then z is not is the residual spectrum of Pi, i — 1,2. 

5.3. Results. 



The orem 5.1. [15] Consider two operators Pi : B — > B which satisfy (5.2), (5.S) 
and §5^. FolQr € (a,M), let 

"1 = h — H 

hxr/a 

AnSC 3 C 2 (C 2 + C s + 2)MH Str (P 1 ) \n(M/r) 

n 2 = ; 1 hni- -. . 

mr/a mr/a 

If 

IIIP1-P2IH < 7, — :=ei(i\,r,£) 

111 2111 -8C 3 (HsAPi)C 3 +jg±- r ) U ^ 



i/ien /or eac/i z e C \ V^ r (Pi), we /law 

iic.-«™is«£±^£)' 



1 



2ei 



Sef 



ln(r/a) 

7 = 



ln(M/a) ' 

8M(C 2 + C 3 ) 2 (M\ ni ro d 2d 



[2C 2 (C 2 + C 3 



M-r \ r J L v ' M -r J M -r' 

& = -rr [MC 2 {C 2 + C 3 + 2) + C 3 ] (C 2 + C 3 ) 2 — + 2C 3 , 

M — r \ r J 

1( ^In this paper will be interested in those r £ («, 1). 
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and 

£ 2 (Pi,r,£) 



-f-V 

tC, \MJ 



4C 3 V M) \ H s . r (P 1 ) [C 2 (C 2 + C 3 + 2)M + C 3 ] + 2C 2 (C 2 + C 3 ) + ^ t 
If 

(5.5) 1 1 \Pi - P 2 | 1 1 < min{ £l (Pi,r, J), e 2 (Pi, r, <5)} := £ (Pi , r, 5) 
i/ien /or eac/i z G C \ Va !r (Pi), we /law 

(5.6) IIKz-P^- 1 -^-^)- 1 !!! < |||P 1 -P 2 ||r(a||(z-P 1 )- 1 || + fe||(z-P 1 )- 1 || 2 ). 

Corollary 5.2. [TS] // 1 1 |Pi -P2 1 1 1 < £i(Pi,r,(5) thena(P 2 ) C Va>(Pi). in addition, 
if \\\P\ ~ P2IH < £ o(Pii r i^)7 i^eri wi eoc/i connected component ofVs tr (Pi) that does 
not contain feoi/i o"(Pi) and cr(P 2 ) have the same multiplicity; i.e., the associated 
spectral projections have the same rank. 

6. Computing the decay rate for £ 
6.1. Applying Keller-Liverani's theorem to C and C m . 

Lemma 6.1. Let denote either C or C m and let f E C Lvp (I). We have 

(1) for all n > 1, |£™| < M n , where M = B Q + 1; 

(2) for all n>l, \\£^f\\ < a n \\f\\ + C 3 M n \f\, where C 3 = + 1. 

(3) |||(£-£ m )||| < Te, where T = 2 • max{a + M, B±}. 

Proof. The inequalities in the first two statements for C are dominated by those of 



C m . Therefore we prove them only for = C m . For |/| < 1, using Lemmas 3.1 



and 4.3 we have 

\c m f\ = |n m £n m /| < < m. 

This implies the first statement of the lemma. To prove the second statement, we 
observe that 

Lip(£ m /) = Lip(n m £II m /) < Lip(£Il m /) 
l<, '' ! <aLip(n m /) + P!|n m /| < aUp(f) + Bt\f\. 

Consequently, 

Lip(CJ) < a n Up(f) + (1 + a + ■ ■ ■ + a n ~ 1 )M n - 1 B l \f\ 
and, for all n > 1, 

\\C n m f\\ <a"||/||+C 3 M«|/|. 



Finally, using Lemma 4.3 and inequality (6.1 1, we obtain 

IOC - c m )f\ < |(£ - n m £)/| + |(n m £ - n m £n m )/| 

< 2eLip(£/) + 2eMLip(/) < 2s[(a + M)Lip(/) + P x |/|] < Te 



a 
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6.2. £ as a perturbation of C m . All the constants arising in Theorem |5.1| are 
(in principle) computable for the matrix representation of the finite-dimensional 
operator C rn , including the number H$ tr (C m ). Thus, using the idea of |17j . we are 



going to apply Theorem 5.1 with C m as P\ and C as the perturbation P%. This 



entails some a priori estimates. 

Lemma 6.2. Given C, 5 > and r G (a, I), there exists £3 > such that for each 
r\ with < mesh(rj) < £3, we have 

(6.2) mesh{ V ) KT^eoiC^r^), 
and 

(6.3) |||£ m -£||| <e (P v ,r,S). 

Proof. See Lemma 4.2 of [TTj^ □ 

In the next lemma we provide a computable upper bound on H$ r (£ m ). 
Lemma 6.3. For f G (C, || • ||) with \\f\\ — 1 and z e C \ Vs,r{C m ) we have: 

(1) ||(z - c m yif\\ < + i)|(z - r™)- 1 ^/! + ^ + f ; 

(2) Kz-r™)- 1 ^,/! < ||(^-£ ro )-i/|| + f. 

Proof. We have 

(z £m) ^ (^ £771) £m ^ 

Then using the fact that IT m is a projection and statements (1) and (2) of Lemma 
|4.1| we obtain 

IK* - c m )-\f - n m /)|| = \\{ z -\z - c m )- l c m + z- 1 ■ i)(f - n m /)|| 
(6 ' 4) ^ll^-^IK A|| fM < 2 . 



z z \z\ r 

Now, write [z— £ m )~ 1 H m f = h, then h — -(£ m h + H m f). Therefore, by inequality 



(6.1 ), we have 



Lip(fc) < r^(Lip(£ m h) + l) < -(aLip(h) + B 1 \\h\\i + l). 
\z\ r 



Hence, 



(6.5) | N |<(_^ + 1)^1 + -1 



and by (6.4) and (6.5) the first part of the lemma follows. For the proof of the 



second statement, by (|6.4), we have 

(6.6) \(z - £ m )" 1 n m /| < - L^-^JW <\\(z- £ m )-7ll + -. 

r 



□ 



^Although the norms used in |17| are BV and L 1 , the proof of his Lemma 4.2 is valid in the 
general setting of 1151 and does not really depend on the norms involved. 



12 Rigorous Approximations In The L°°-Norm 

Now we define 

H* Sr (C m ) :=sup{(^ Bl ^+l)|(z-£ m )- 1 n m /|+^^ + 2 : ||v|| x = 1, z e C\V s , r (£ m )}, 
r — a r — a r 

t ln8C 3 C 2 (C 2 + C 3 + 2)MH* Sr (C m ) \n(M/r) 1 

n 2 '■= i / ■ \-ni- — 

In r/a mr/a 



8C 3 (Hg r (£ m )C 3 



M-r> 



sl[C m ,r, 5) 

L f-V 

Co \MJ 



AC 3 KM J \Hl r (£ m )[C 2 (C 2 + C 3 + 2)M + C 3 ] + 2C 2 (C 2 + C 3 ) + ^ 
and 

e* {C m ,8,r) := min{el{C m ,r,5),e* 2 {C m ,r,5)} 

Note that H^ r {C m ) provides a computable upper bound on Hg tr (C m ), and conse- 
quently £g(£ m , 5, r) provides a computable lower bound on eo(C m , S, r). The critical 
step is to obtain control on the separation of the point spectrum of C outside the 
essential spectral radius a. More precisely, the following algorithm will compute 
numbers 8 — 5 C , a < r = r c < 1 such that S c < 1 — r c and e = e c > such that for 
any r\ with mesh(r/) = e c 

(1) meshfa) < (r^eotA,,,^); 

(2) B(p, 5 C ) n B(0, r c ) = 0, where p is the dominant eigenvalue of C m . 

(3) If pi 7^ p is an eigenvalue of C m , then pi 6 P(0, r c ). 

Thus we obtain the a spectral gap for C m and consequently, by Theorem |5.1| a 
spectral gap for C. 

Algorithm 6.4. Given T and p€N, then 

(1) Picfc 5=1/ j> 

(2) 5e* r = 1 - 25. 

(3) Choose m, the number of partition intervals. 

(4) Set e = i £/ie mes/i size 0/ i/ie partition. 

(5) find i/ie matrix representation of £ m 

(6) C/iecfc i/e < (r)" 1 ^^™, <5, r). 

If (6) is not satisfied, feed in a larger m repeat (3)- (6); otherwise, continue. 

(7) Check that B(p,5) n B(0, r) = and pi £ P(0, r) /or any eigenvalue pi of 
C m with Pi ^ p. 

(8) // is satisfied, report 5 C := 5, r c := r and m c := m; otherwise, multiply 
p by 2 and repeat steps (l)-(l) starting with the last m that satisfied (7). 

Proposition 6.5. Algorithm\6.J\ stops after finitely many steps. 



Proof. The proof is similar to the proof of Proposition 1 of pQ 



□ 
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6.3. A computable bound for the rate of decay of correlations of L. We 

will benefit from Algorithm |6.4| in many directions. First it provides us with a nice 
spectral picture of £; i.e., an estimate on the size of its spectral gap; moreover 
it enables us to compute an upper bound on the norm of the resolvent of the 
continuous (infinite dimensional) operator C bounded away from its spectrum. Such 
a computable estimate was impossible to do before Algorithm |6.4| 

Lemma 6.6. We have 

sup{||(z - : z e C \ V Sc , rc (C m J} < 4(g2 + C3)(M/r e )"I ;= 

(M - r c )2e\ 



wh 



ere 



(Tc/M) 

El — 



and 



n , = ln2C 2 
1 lnr c /a 

r m8C 3 C 2 (C 2 + C 3 + 2)MH* r (£ m J ln(M/r c ) n 

1 — / ~ + n i n — / — I ■ 

mr c /a mr c /a 



Proof. We have e c < T 1 eo(£ mc ,S c ,r c ). Therefore, by Theorem 5.1 

K5 c ,r c (£m c ), 

< 4(C2 + C3)(M/r)"i < 4(C2 + C3)(M/r; 



»i 



(M-r c )2e 1 (£„ ic ,(5 c ,r c ) " (M - r c )2ej 

Here we have used the fact that Hs c ^ c (C mc ) < Hg r (£ mc ) which implies that 
Ei(£ m< ,, S c , r c ) > e|. This completes the proof of the lemma. □ 

We now provide a computable estimate on the rate of decay of correlations of C. 

Lemma 6.7. Set R := min{(l-5 c +r c )/2, (p c - 5 c + r c )/2} . Then for f e C Llp (I), 
J f = 0, we have 

\\C n f\\<R n+1 Hl rr (C)\\f\\. 
Proof. Let t^r be the following spectral projection 

7TR := -^r / (z- Cy 1 dz. 

lm J{zeC: \z\=R} 

Then for any / € C Up (I) 

Cf = f* I f + Cn R f. 



This implies that if J f = we have 7Tr/ = / and consequently 

1 
2iri 



C n f=^-[ z n {z-C)- l dz. 

l{z£C: \z\=R} 



Thus, 

\\£ n f\\<R n+1 sup \\(z-C) 

\z\ = R 



< sup IK*-/:)- 1 !! Il/ll < R n+1 H^ rc (c)\\f\\ 

{zee\Vs c ,r c (£ mc )} 



a 
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7. Rigorous approximation with fast rates 

The projection-based discretization does not appear to be particularly well suited 
to estimating the rate of approximation of the Perron eigenvector /* m for C m to 
the invariant density /* for C, where rj m denotes the uniform m— cell partition. 
We discuss this in a little more detail in the next section. However, we can use 
the spectral estimates on C obtained via perturbations C rn combined with the 
discretization scheme associated to P m to obtain good, rigorous approximation 
rates in the uniform norm. 

We begin by deriving a suitable Lasota-Yorke inequality for the operator P m . 

Lemma 7.1. For f e C Llp {I) we have P m / G C Lip (I) and 

Lip(T m f) < pLip{f) + 
where [3 = (3/2)a and B 1 = (Z/2)B X . 

Proof. The proof follows from Lemma |3.3| and |4.4| □ 
Without loss of generality, we will assum^] (3 < 1. 

LelQ f* m denote the nonnegative, normalized Perron eigenfunction of C m corre- 
sponding to the eigenvalue 1. Our goal is to prove \ff n — f*\ = 0{ Xj ^) 

First, we find a computable upper bound on the Lipschitz norm of f* n . Due to 
the fact that \f* n \ > 1, this is best done in two steps. 

Lemma 7.2. and f* satisfiy the following estimates 

(1) \fa\<Ki, 

(2) |/*| <K 2 , 

(3) Lip(f* m ) < j^Kl 

(4) Lip(f*) < ^K 2 

where K x := ^ + 1 and K 2 := (^) + 1. 

Proof. For the first estimate, note that the matrix representing P m (as determined 



in Section 4.2) is stochastic and strictly positive, hence mixing. It follows that 

f* = limP n 1 

J TIT, LULL 1L m _L 



uniformly on the interval /. But, as a consequence of Lemma 7.1 we get 

Lip(P^l)| < ^Lip(l) + ^111 = Y^-p- 

Since P"„l is nonnegative, with unit L 1 — norm we find |P"jl| < 1 + jzV This is 
the first estimate. 

The second estimate is similar, using Lemma |3.1| in place of Lemma |7.1| to 
obtain a uniform bound on the Lipschitz constants and L°° — norm of the iterates: 
Lip(£"l) < an d |£™1| < ^ + 1. By Helley's Theorem, we may obtain a 
uniformly convergent subsequence from the sequence of iterates; the averages along 
this subsequence converge uniformly to /*, giving the second estimate. 

Next, using Lemma|7.1| we obtain 



Lip(/* ) - Lip(P m /* ) < pUp(f^) + Bil/* |. 

12 In fact we can choose kg > k so that that MaJ < 2/3 which makes < 1. 

^Note the difference between and /* m the Perron eigenvectors for P m and Cm respectively. 
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Therefore, 



The same argument, using Lemma |3.l| and estimate (2) gives the inequality in part 
(4). ' ' □ 



Up(/*)<^|&|<^*, 



Now we state and proof of our main result, Theorem |7.3| There are certainly 
many ways to proceed which lead to a convergence rate of order ln ^"^ p*j However, 
different proofs will lead to different constants multiplied by this rate. Since our 
target is to produce explicit computation of the error bound, we write the proof in 
a way which makes the constants as small as possible. 

Theorem 7.3. For any m £ N, we have 

{Kl + Kl)RHl rc (C) + ^rin(m)/ In iT 1 !^! 



\r-a< 

Proof. We have 

\r - a < rx - £ n a + \£ n r m - £ n r i 
( ' } =(i) + (ii). 

Using Lemmas |6.7| and |7.2| for any n we obtain 

(7.2) (II) < (7^(1 + + K 2 (l + ^))Hi^{C)R n+l . 

In (J) the problem is that we only have the weak estimate |PJ^| < (2M)™. However, 
we know from Lemma 3.1 that for all n > 1, \C n \ < M. Therefore we change the 
order in (I) and benefit from Lemmas |3 . 1 1 and 7.2 So, we write 

'9-1 f * I 
mj* m J ml 



(I) = \C n f m - ¥ n m f m \ < J2 \£ n ~ q (£ - Pm)P^ 1 /, 



q=l 

n n 

(7.3) < m£|(£ - P m )(P?- 1 /;„)l = m£|(7 - Q m )(£P«T 1 /* 

q=l q=l 

5M 2 A„, 5M 2 T ^ 

q=l 



Choosing n = [ln(m)/ln(i? 1 )] and using (7.2) and (7.3 1 completes the proof. □ 

8. Example 
In this example we use the map: 

(8.1) r(x) -- 



£f for < x < ^ 



12a; - i for < x < ^ ' 



where i = 1, .... 11. 



To improve the rate of convergence from Q( lnm ) to O(A), one would need Lip(/* — 
Qmf*) = 0( — ). However, this cannot be obtained unless /* is very regular. For the class of maps 
under consideration /* is, in general, only Lipschitz continuous and hence Lip(/* — Qmf) = O(l). 
The situation is analogous to the failure to obtain rate O(i) for Ulam's method via spectral 
perturbation arguments; see \7\. 
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8.1. Computing Hg r (£). We will use Equation 3.1 to estimate the relevant 
constants. For / G BV we have 

V r £ T /<l/liy/ + 2/ll||/|| 1 . 

Consequently, for any n > 1, we obtain 

V£?f< (l/ll)"V7 + l/5||/||i. 

Hence a = 1/11, B = 1/5, and M = 6/5. We set T := r. Then a = a M = 6/55. 
For / G C Up {I) and any n > 1 

Lip(£/)<6/55Lip(/) + 2/(ll)Vl- 

We now use Algorithm |6.4| to obtain rn c which allows us to compute an upper 
bound on the norm of the resolvent of the continuous operator Hg r (CjQ 



Table 1. The output of Algorithm 6.4 for the map in (|8.1 



r 


0.8 


0.8 


6 


0.1 


0.1 


e 


2 x 10" 4 


8 x 10" 5 


H* Sr (C m ) 


27.2822974 


27.29122985 




0.0001276212886 


0.0001275730002 


Loop I 


Fail: reduce e 


Pass 


Loop II 




Pass 


Output 




mesh(m c ) := 8 x 10" 5 , 5 C := 0.1, r c := 0.8 



Using this output, by Lemma |6.6| we obtain 



H, 



(£) := 



4(C*2 + C*3)(M/r c )™ 1 



< 36366.11326. 



(M - r c )2e\ 

8.2. Computing a uniform bound on the approximation error. The first 
few eigenvalues of C are estimated (to 4 decimal places) by 

1.0000, 0.0901, 0.0088 0.0000 + 0.0069i, 0.0000 - 0.0069i 



Next we compute the constants from Lemma 6.7 and Lemma 7.2 

521 549 

K 1 = ,K 2 = and R = 0.8500 

506' 539 

Then in Table [8T2| we list the number of bins m and the corresponding approximation 
error of \f 



The results in Table 8.2 are undoubtebly far from optimal. We 
Table 2. Error := \f* 



fm\, vs m where /* is the invariant 



density of the map in (8.1 ) and is the invariant density of £ r . 



m 


Error 


10 y 


0.06515863735 


7 x 10 6 


0.009314201609 



reiterate that efficiency in computations has not been our main focus. 



^For details on how to use the computer to compute an upper bound on Hg r {Cm) in the 
algorithm for the discrete operator C m see pQ. 
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9. Discussion on the rate of convergence of the discretization 

SCHEMES C m 

The operator C rn is not Markov and its dominant eigenvalue p is typically > 1. 



This hinders an estimate similar to that of (/) in Theorem 7.3 In particular the 
sum J2q = i£-mfri m cann °t be controlled since £ m f* m = p q f* lm - Thus the key trick 
for obtaining rate of convergence elne -1 is unavailable for the scheme H m (unless 
it happens that p < 1). 

In principle, one can obtain rigorous, computational estimates for the scheme 
Tl m , through direct application of Theorem |5.1| obtaining L°°-norm estimates in 
terms of the difference of the spectral projections 



iProji - ProjJ <5 c [aH* 



bH* 2 \e 



7c 



\n(r c /a) 

\u{M/a) 

to achieve an error smaller than 10 -2 , one would need e < e 



< 1. From this estimate, it follows that for the map in (8.1|, 

35.9 



where 7 C = 

Needless to say, 
this is not a practical approach in general. 

On the other hand, one expects that the both discretization schemes should 
be more efficient than indicated by the above theoretical estimates. We conclude 
this article with a simple numerical example (non-rigorous) that shows the kind of 
performance one should be looking for, even in the non-Markov discretization case. 

dens_blue approx_plus_red 




Figure 1. The graphs of /* (blue) and ft + 0.02 
(red) for II m — scheme 
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Consider the following Markov interval map which is only eventually expanding: 



(9.1) t(x) = \ \ 



2a; £„_ n ™ 1 



2.r 



for < x < 3 
for \ < x < 1 



and whose invariant density is given by 

Since the invariant density /* is known for this example, we can check the efficiency 
of our discretization scheme U m . Using m = 10, we plotted in Figures [T] and [2] the 
graphs of /* , /* +0.02 and f* — f^ m respectively. As the figures show, our tolerance 
A is attained with only 10 bins; the efficiency of the discretization scheme II m in 
achieving an L°° approximation is pretty impressive. 




Figure 2. The graph of /* - 
for the II m — scheme 
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